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Within the frame of a novel treatment we make a complete mathematical 
t> \ analysis of exactly solvable one- dimensional quantum systems with non-constant 



mass, involving their ordering ambiguities. This work extends the results recently 
reported in the literature and clarifies the relation between physically acceptable 
effective mass Hamiltonians. 
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The study of quantum mechanical systems with position dependent mass has 
been raised some important conceptional questions, such as the ordering am- 
biguity of the momentum and mass operators in the kinetic energy term, the 
q-i boundary conditions at abrupt interfaces characterized by discontinuities in the 

mass function, etc. Therefore, the form of the effective mass Hamiltonian has 
been a controversial subject in the literature. In recent years there has been a 
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growing interest in the study of such systems due to the applications in condensed 



matter physics and other areas involving quantum many body problem. These 
applications stimulated a lot of work in the literature regarding the development 
of techniques for the treatment of such systems, for a recent review, see [1-6] and 
the related references therein. In all these works the main concern is in obtaining 
the energy spectra and/or wave functions for quantum systems with spatially 
dependent effective mass. Moreover, exact solvability requirements result in con- 
straints on the potential functions for the given mass distribution. Though there 
has been a large consensus in favor of BenDaniel and Duke Hamiltonian (BDD) 
[7] proposed in the literature as an appropriate one, the question of the exact 
form of the kinetic energy operator is still an open problem for such systems. 

Within this context, the present Letter involves an alternative scheme to ob- 
tain unambiguously the Schrodinger equation with non-constant particle mass, 
which makes clear the relationship between the exact solvability of the Schrodinger 
equation and the ordering ambiguity. The model explored here restricts naturally 
the possible choices of ordering and provide us a clear comparison between the 
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solutions of different but physically plausible effective Hamiltonians clarifying the 
physics behind ambiguity. 

To achieve our goal defined above, the recently developed non-perturbative 
technique [8] is employed within the frame of supersymmetric quantum mechanics 
[9]. In this unified model, the BDD Hamiltonian is considered as an unperturbed 
term while modifications due to other effective Hamiltonians are treated as an 
additional potential in the same framework. This realization is of prime sig- 
nificance in the calculation of physical processes, which so far did not receive 
adequate attention. 

In the following sections the model used through the work is introduced and 
its applications are presented where the superiority of the present scheme is also 
discussed. 

There are several ways to define the kinetic energy operator when the mass 
is variable. Since the momentum and mass operators no longer commute, the 
generalization of the Hamiltonian is not trivial and this kind of physical problem 
is intrinsically ambiguous. Starting with the von Roos effective mass kinetic 
energy operator [10], which has the advantage of an inbuilt Hermicity, 

H uR = - [m a {z)pmP (z)pm' 1 (z) + rn> \z)finft \z)pm a (z)] + V(z), (1) 

where a+P+j = —1. By the correspondence in wave mechanics p — > —ih-f^, z — > 
z and on setting 

m(z)=m M(z), % = 2m Q = l, (2) 

where M(z) is the dimensionless form of the mass function, the effective mass 
equation can be written in a differential form, 



d 

dz 



1 <iNf(z) 



M W dz j +»""«*<*> = **<*), (3) 

Here, V e ^(z) is termed the effective potential energy whose algebraic form de- 
pends on the Hamiltonian employed 

V<"(z) = V (z) + U ai {z) = V (z) - (^)^ + (ary + a + , (4) 

in which the first and second derivatives of M(z) with respect to z are denoted by 
M' and M" , respectively. The effective potential is the sum of the real potential 
profile Vq(z) and the modification U ai (z) emerged from the location dependence 
of the effective mass. A different Hamiltonian leads to a different modification 
term. Some of them are the ones of BDD (a — 7 = 0), Bastard [11] (a = — 1), 
Zhu-Kroemer (ZK) [12] (a = 7 = -\) and Li-Kuhn [13] (/3 = 7 = -\). 

Considering the supersymmetric treatment of effective mass Hamiltonians by 
Plastino and his co-workers [14] 

1 d%> d ( * \ 

Am = -^^- + W^, A + * = — +WV, (5) 
y/M dz dz \Vm) 
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where A and A + are linear operators and W(z) is a superpotential, the super- 
symmetric Hamiltonians are expressed as 



Hl = A+A= i +w *.(w), (6) 

Mdz 2 \MJ dz \VmJ 

and 

From which, supersymmetric partner potentials are 

V susy = W 2_(W\' SUSY = SUSY + 2^_[ 1 )( 1 V_ (g) 



At this stage, we use the spirit of recently developed non-perturbative approach 
[8] by expressing the total wave function as a product, 

q?( z ) = $(z)G(z). (9) 

In the above equation, $ denotes the wave function corresponding to the un- 
perturbed piece of the effective potential in Eq. (4) while 6 is the moderating 
function due to the modified term U aj therein. 
The use of (9) in (3) yields 

T, (t + e" + 2 ¥ e J - Jfl { ¥ + e J = " (10) 

which reduces to the usual Schrodinger equation with a constant mass when 
M — > 1. With the consideration of (6), where the superpotential now can be 
given as 

W(z) = W (z)+AW(z), (11) 

with Wq and AW being superpotentials corresponding to the unperturbed po- 
tential (Vo) and modification term (U aj ) respectively, Eq. (10) is transformed 
into a couple of equation, 

^-(Si^-** w ' = -jnT' (12) 



Aly2 - (%)' + 2W ° AW = u " - AE ' AW = -^1- (13) 



In the above equations, E = E + AE due to V e ff — Vo + U aj . Therefore one can 
easily see the contributions, if any, to the energy and wave function due to the 
use of effective Hamiltonians other than BDD which represents the unperturbed 
Hamiltonian in the present scenario since it has no modification term, see (4). 
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We are familiar with (12) as a standard supersymmetric treatment of the 
Schrodinger equation for the exact solutions. However, Eq. (13) is new and 
is the most significant piece of the work presented in this letter. Because it is 
a non-perturbative approach by Riccati equation, which reproduces the whole 
corrections coming from U aj if, of course, Eq. (13) is exactly solvable. 

To proceed we remind a general consensus [5, 13] that the resolution of the 
ordering ambiguity in this problem could come from a scheme that starts with 
the relativistic Dirac equation with spatially varying mass then taking the non- 
relativistic limit. This is due to the fact that the Dirac equation is inherently 
free from the ordering ambiguity and that taking the non-relativistic limit is a 
well defined procedure. Bearing in mind this point we propose a correct choice 
of AW as 

which directs us to find correct ordering parameter(s) leading to the physically 
plausible effective Hamiltonian(s). Through Eq. (13), the parameters get decou- 
pled in a natural way and the ambiguity in the choice of proper kinetic energy 
operator disappears. Substituting (14) into (13), we obtain 

( AW\' 

AW — i-y=\ = U aj , AE = -2W AW, (15) 



if either a = 7 = which yields the BDD Hamiltonian or a = 7 = — \ corre- 
sponding to the ZK Hamiltonian. It is stressed that the results are independent 
of any choice of M(z) and in case a = 7 = Eq. (13) vanishes. This restriction 
is in agreement with the discussion in Ref. [15] and also with the work of Bagchi 
et all [3]. 

Though the present formalism has a wide spread applicability, for clarity we 
now simply consider the two examples which were investigated in Ref. [14]. This 
consideration will shed a light in understanding the interrelation between the 
BDD and ZK effective Hamiltonians bearing in mind the results presented in [14] 
for the systems of interest. 

The simplest case of the shape invariance integrability condition [9], leading to 
exactly solvable potentials, corresponds a uniform energy shift e between partner 
potentials, 

V 2 SUSY (z, e) - Vf^iz, e)=e = 2E (16) 

since AE term appearing in the partners due to U aj cancels each other. The 
replacement of (8)into (16) gives 

2 W+ Avn ( 1 \ ( 1 V' =Ei (17) 



M 
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from which one finds the superpotentials leading to the hamiltonian with Vq, 



since AW = -(M'/2M 3 / 2 ) = (lfy/M^. To finalize the full treatment, one needs 
the total superpotential, W = Wo + AW from which the results in [14], Eq. (35) 
and the subsequent equations, can easily be reproduced. 

From this short discussion, it is obvious that (i) there will be no contribution 
to E Q due to the modification term. For this reason total energies in both system 
having a constant mass and position dependent mass are equal, (ii) From (13), 
the contribution of U aj to the unperturbed wave function is (for the ground state) 

Q n =o(z) = exp (- f ^M(y)AW{y)dy} = m 1 ' 2 . (19) 

Thus, going back to (9) along with Eqs. (12) and (18), the full unnormalized 
ground state wave function is expressed as 

$n=o(z) = \m~ l/ \z)§(z)\ m 1/2 (z) = m l/ \z)§{z), (20) 



where z = J z y M(y)dy, which supports the reliability of the present formalism 
[1]. The excited state wave functions can be determined [9] in algebraic fashion 
by successive application of the linear operators in (5) upon the ground state 
wave function, (iii) The both choice, namely the BDD and ZK Hamiltonians 
are represented with a unique superpotential leading to exactly equivalent wave 
functions for shape invariant potentials, (iv) From (8), as a — 7 = — \ , one gets 

2W' 



ySVSY = (ySVSY + ^ + ^_ ^ 



pointing a duality between BDD and ZK schemes, which reveals the suggestions 
in [1,3]. 

Let us proceed with another example in Ref. [14] where the superpotential 
leads to a Morse-like spectra, 

W(z,A)=A + f(z), (22) 

in which, within the frame of the present formalism, f\z) = }'q{z) + Af(z) that 
turns the form of (22) into 

W(z, A) = [A + f (z)} + Af(z) =W + AW(z) (23) 

From the shape invariance condition V2 USY (z, A) = Vf USY (z, A-\)+R(A) used 
in the supersymmetric quantum theory [9], where A is the potential parameter 
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and R involving both parameter, A and A , leads to the ground state energy of 
the system. In the light of the work carried out in [14], the substitution of (23) 
in (8) within the frame of shape invariance condition above produces 



2 Wo + A/') 



M 



M \VM 



= A 



M , 



2A (/o + Af) . 



(24) 



Remembering AW = Af = (-^7=g) for a = 7 = — |, the above equation 
rearranged as 

f^z) + b 1 (z)f (z)=b 2 (z), 



IS 



where 



61 = \VM, b 2 



A 



M 



M , 



+ 



M, 



(25) 

(26) 
-b 2 . 

(27) 



where C is an integration constant. Employing the mass function used in [14], 
M =[(a + z 2 )/{\ + z 2 )} 2 , we obtain 



From (24) it is clear that Af term affects only b 2) since when Af — > b 2 
The solution of differential equation in (25) gives 

f Q (z) = \c + J b 2 (y)dy exp J &i(f)df J x exp - J bi(y)dy 



W(z) = W + AW 



( z(a — 1) \ z(a — 1) 

= [A + Cexp [-A {z + (a - 1) arctanx}] - -± -L + 2-^ -L, (28) 

V (a + z 2 ) ) (a + z 2 ) 

that is Eq. (53) in [14]. From (12), the corresponding potential function, energy 
and wave function can be expressed as in [14] , which are out of interest in this 
letter. Generalization of the above discussion to a formalism which is applicable 
to all spatially varying masses, yields 



W(z) = Wq+AW = {A + Cexp 



-J h(y)dy 




Plastino and co-workers [14] studied this problem in case a = 7 = considering 
only the BDD Hamiltonian and arrived at Eq. (53) in their work, which addresses 
(28) in our work. This means that BDD and ZK effective Hamiltonians in fact 
reproduce same results employing an identical superpotential, which once more 
supports the realization introduced by (21) that they are their supersymmetric 
partners. 

In this work we have discussed the problem of solvability and ordering am- 
biguity in quantum mechanics for the systems with a position dependent mass. 
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The present scheme restricts the possible choices of ordering. Proceeding with 
this consideration it has been observed that the only physically allowable BDD 
and ZK Hamiltonians are in fact their supersymmetric partners that reproduce 
identical results in their independent considerations due to use of an identical 
superpotential. We hope that this observation would make a contribution to the 
ongoing debate in the literature regarding the isospectral effective mass Hamil- 
tonians. 

The authors wish to thank the referee for his helpful comments and sugges- 
tions. 
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